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Abstract 

In this work we solve the Dirac equation by constructing the exact 
bound state solutions for a mixing of generalized vector and scalar 
Hartmann potentials. This is done provided the vector and scalar 
potentials hold some relation. Namely, one must be equals to or minus 
the other. Finally the case of some quasi-exactly solvable potentials 
are briefly commented. 
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Very recently in this journal C. Y. Chen [T] presented exact solutions for 
fermions in the presence of a classical background which is a mixing of the 
time-component of a gauge potential and a scalar potential. He has con- 
sidered the particular case where both potentials are of the Hartmann type 
potential [7]. In fact, he in collaboration with other authors has dedicated 
great attention to this mater in the recent years, which can be verified from 
a number of interesting works On the other hand, many years ago, 

Hautot [8] solved the Schrodinger equation for the Coulomb and harmonic 
oscillator potentials with some asymmetrical terms like added. Here 
we intend to show that one can find some other cases for which the Dirac 
equation with classical potentials of vector and scalar natures with spherical 
asymmetry can be solved exactly. We begin our study treating the case of 
the Hartmann-type and after that we discuss a generalization of Morse-like 
potentials [3]-[T0j with spherical asymmetry. Finally we comment on the 
quasi-exactly solvable potentials p!T]-[T2]. 

The generalized Hartmann potential is defined here as 

The time-independent Dirac equation for arbitrary scalar and vector poten- 
tials looks like 

ca-P+l3 {Mc^ + S{r)) 
where it is defined that 

with a the vector Pauli spin matrix and / the identity matrix. Now, using 
the Pauli-Dirac representation, with 

(4) 

we get the following set of coupled equations for the spinor components 
ca-Px{r} = [E-V{r}-Mc^-S{r}]^{f), 

(5) 

ca-Pip{r) = [E-V{r) + Mc'^ + S{r)]xir)- 



^Ij{^ = [E-V (f)] ^ (f) , (2) 
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At this point we can treat two non equivalent exact situations. The first 
when S {'f) = V (r) , which was the one considered in [T] , and another where 
S {-r) = —V (r). Once the treatment is quite similar in both cases, we start 
by dealing with the first one and then comment about the second case. 

The case with S {r) = V (r), allow us to decouple the Dirac equation as 



ca-P 

E + Mc^ 



c^P 2 + 2 {E + M c^) V (r) ^ (rj 
leading us to the following Schrodinger like equation 

'v,\ fie) 



[E^-M^c^]^i 



(6) 



(7) 



-fi'c^V^- {E + Mc^) 



^ (f) = (E' - M'c^) ^ (f) . (8) 



Performing now the usual separation of variables in spherical coordinates 

if (r) = ^= ^ e (0) , m e (9) 



^2ti r 

we obtain the equations for u (r) and G [9) , which are respectively given by 
1 



sin {9) d9 
and 

d^u (r) 
dr'^ 



d f ^d&(9) 



m 



sin' 9 



+ {E + Mc')f{9)-s 



e{9) = 0, (10) 



+ 



E + Mc^ VqX 



u (r) 



u (r) 



From the above equation one can see that bound-state solutions are pos- 
sible only ii \ E \ < Mc^ and that there is no room for bound-states for 
s < —1/4, because this is the critical value of this parameter, due to the fact 
that below this value, the singularity of the potential produces the so called 
fall to the center. Finally, regarding the product of parameters Vq A, there is 
only one situation where there are no bound states. That happens if Vq A < 
and s > simultaneously. 
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Now, we are able to specify the form of the function f{6), in order to 
guarantee that fllOp can be exactly solvable. Following the original work of 
Hautot fU, we can consider the cases where 



(7 + /3 cos 9 + a cos^ 
sin^ 9 

(7 + /5 cos^ 9 + a cos^ ' 



/iW = ^^^^ ,V ^ (12) 



= . 2n 2n (13) 

sm 9 cos^ 9 

f^{9) = j + p cot9 + a cot^9. (14) 

For each one of the above configurations for f{9), Hautot was able to 
map equation ffTOl) into a hypergeometric differential equation, which has 
finite solutions for G {9) in the range < < vr provided that s satisfies 
some constraint as a function of the parameters a, f3, 7 and a new quantum 
number k = 0, 1, 2, that was introduced in order to render (9) a finite 
polynomial (Jacobi polynomial). Here we present a detailed analysis of the 
first case and comment about the other cases. 

For fi{9) and f2{9) the regular solutions found for 6(6') are generally 
given by 

Qk{0) = z^{l-zY,F2i-k,b,d;z), (15) 

where iF2{—k, b, d; z) is the hypergeometric function and z = cos^ 9 
For fi{9) one has 

P = l[m' + {E + Mc'){a-P + j)Y^\ 

V = l[m^ + {E + Mc^){a + P + ^)Y^\ (16) 



b = k + 2{p + u) + l, 

d = l + 2p, (17) 

and 

s + \ = \ib + kf - {E + Mc^)a. (18) 

At this point some discussion on the range of validity of the potential 
parameters should be done. Note that in order to keep b and d real, it is 
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necessary to impose that a+'j—sign {f3) , (3 >0. Furthermore, the imposition 
that the system should avoids the fall to the center, implies into the following 
additional restriction. 



+ 



1 



{E + Mc^)a > 



16' 



(19) 



Defining now the variables y = VE + Mc^ , 5i 
q; + /3 + 7, we obtain the following equation 



a 



/? + 7 and ^2 



(\/^+^)'-4a |/' + 2 (^v^+ v^) 1/ + ^ >0. 



(20) 



On the other hand, we know that y is a. positive definite variable. As a 
consequence, it is not hard to conclude that the only way to avoid any further 
restrictions over the energy of the system, is to require that the coefficient 
of be positive definite also. This requirement implies into the following 
equation for the parameters 



a + 7 + \/(« + 7)^ + /^^ - 4 a > 0. 



(21) 



A very similar analysis could be done for /2(^), because in this case we 
have 



\ + \[l + A{E + M^)^]"\ 



(22) 



h = /c + 2(p + i/) + ^, 

d = 2p+^, (23) 

and 

s + 1 = (6 + A;)2 - (^ + M^)a. (24) 

Its analysis is straightforward and we let it for the interested reader. 
Althought the case with needs a different change of variables in order 



5 



to cast the equation (fTUj) into a familiar hypergeometric one, it presents the 
same general behavior in terms of restriction over the potential parameters, 
and we do not present it here also. 

It is worth to mention that our analysis could be extended by including 
some other cases for f{6), as was done in reference [8j, had we taken into 
account a two-dimensional scenario. 

The solution for the radial equation (|TT1) . by its turn, can be obtained 
exactly following the reference [13]. It corresponds to an effective radial 
equation when a non-relativistic particle of effective mass Me// = 1/2 is 
under the action of an effective potential 

hh E + Mc'VoX 
^-ff = 7^ -2 — (25) 

and effective energy given by 

^^ff -2 ■ (26) 

The normalized radial eigenfunctions u{r) can also be read off directly 
from reference [13j and is given by 



Un,i{r) 



:2k) 



T(n + 1] 



2n[r(n + 2/ + 2)]^ 



N 1/2 

^ exp{-Kr) (2/tr)'+^ L2'+/_^(2/tr) 



(27) 

where L2'+/_^(2Kr) are the Laguerre polynomials, n = 0, 1, 2, ... denotes the 
number of nodes of the radial function. We have defined 



^ ^ and l = 2^/s + l/4-l/2>0. (28) 



The effective energy satisfies the following relation 



{E + Mc^)Vo\ 
2hc{n + 1 + 1) 



(29) 



from which we can take the energy eigenvalues. At this point we could be 
tempted to write a solution 



E = Mc^ 



l + r2 



(30) 



where r 



2hc{n+i+i)' as was done in However, we should remind that 
in this case / is a nontrivial function of the energy leading to a somewhat 
intricate equation for the energy. 



Finally one can note that the equations for the bispinors ip{r) and x(0 
in the situation where S (r) = —V (r) are given respectively by 



ca-P 
E-Mc^ 



(31) 



and 



c2p2^2 {E~Mc^)V 



[E' - M'c'] X I 



These can also be obtained by performing the transformations 



(32) 



V^X, X^^. V^V and E ^ -E (33) 

in the equations (jnD and ([7]). Then the solutions for the situation where 
S {r) = —V {r) can be obtained from those in which S [r) = V (r) by means 
of the set of transformations (!33|) . 
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